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Abstract 

The uncorrelated ("sequential") and correlated ("nonsequential") double ion- 
ization of the H2 molecule in strong laser pulses is investigated using the tools 
of nonlinear dynamics. We focus on the phase-space dynamics of this system, 
specifically by finding the dynamical structures that regulate these ionization 
processes. The emerging picture complements the recollision scenario by clarify- 
ing the distinct roles played by the recolliding and core electrons. Our analysis 
leads to verifiable predictions of the intensities where qualitative changes in 
ionization occur. We also show how these findings depend on the internuclear 
distance. 
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1. Introduction 

The interaction of strong fields with molecules and their molecular ion and 
isoicpic species have been studied extensive* (see Re,s. QQEU reviews). 
Among molecular species, the HVlaser interaction remains a subject of particu- 
lar interest because of the attosecond and femtosecond nuclear time scales and 
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the relatively simple decay channels, showing the way to controlling the motion 



of electrons in molecules 



. Nonsequential double ionization of H2 has 



been studied in detail with an emphasis on the rescattering dynamics (e.g. Q). 

The strong enhancement of nonsequential (correlated) double ionization over 
its sequential (uncorrelated) counterpart is one of the most striking surprises in 
laser-matter interactions. Indeed, at some intensities, double ionization rates 
can be several orders of magnitude higher than the uncorrelated sequential 
mechanism leads one to believe. The hallmark of this strong correlation is 
the so-called knee of the double ionization probability versus the laser intensity. 
An example of a knee is provided in Fig. [1] for a model of the H2 molecule. 
Similar knees have been observed in experimental data §[ {J nj 111 12, 3, 
I4I 0, Q and successfully reproduced by quantal computations on atoms and 



molecules [171 
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. Different scenarios have been proposed to explain 



the mechanism behind this enhancement [8, 
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30|, |31|, |32J, |3 
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recollision scenario 



3 r > 
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0,0, 00, HE, 



361 ] - When confronted with experiments 




231 ]. in which an ionized electron is hurled back at the 



core and ionizes the second electron, seems in best accord with observations. 
It is generally agreed that the sequential process stems from the action of the 
laser pulse, whereas the electron-electron correlation is at the essence of the 
nonsequential process. 

While classical mechanics would, at first sight, seem an inappropriate tool 



for studying this strong-field process, recent classical models 
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40] have been surprisingly successful in capturing key features of 

I I 

the process because of the dominant role of correlation [311 ]. Indeed, entirely 
classical interactions turn out to be adequate to generate the strong two-electron 
correlation needed for double ionization. These ideas have been put to the test 

0, 0, 0, 40 1 . A purely classical scenario 



on the helium atom 
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for the double ionization of the Helium atom has been proposed in Ref. [4C 
based on the identification of the organizing centers of the classical dynamics 
of an effective Hamiltonian model. In the current manuscript, we explore the 
generality of this scenario and in particular how these ideas apply to the double 
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ionization of the H2 molecule driven by a linearly polarized strong laser pulse. 
We shall show that the ideas developed fore helium are also applicable here, 
with some natural variations for the presence of the separated charge centers. 
In summary, the emerging classical scenario for nonsequential double ionization 
is: 

• Without the field, the dynamics of H2 is very chaotic and mainly driven 
by four weakly hyperbolic periodic orbits. These periodic orbits define an 
"inner" and an "outer" electron. 

• At time t = 0, the laser pulse is switched on and focuses its action on the 
outer electron by pulling it out of the core region. 

• After half a period of the field, the outer electron is hurled back to the 
core region and interacts with the inner electron for a very short time. 
The interaction, defined as a recollision, transfers energy from the outer 
to the inner electron. Several outcomes are possible and among them are : 
A double ionization where both electrons leave the inner region, a single 
ionization where one of the electron leaves the nucleus and the other one 
stays in its neighborhood, an exchange of roles between the inner and the 
outer electron, or the outer electron proceeding for another recollision. 

• Between two recollisions, the inner and outer electrons have their own 
respective dynamics. The inner region is mainly influenced by a periodic 
orbit which has the same period as the laser field. As long as this periodic 
orbit exists, it organizes an inner region where all the trajectories of the 
inner electron remain bounded. 

This scenario complements the recollision scenario by adding the phase space 
picture of the inner electron to the one of the outer electron. The mechanism 
in phase space allows us to formulate some verifiable predictions of two char- 
acteristic intensities of the knee : The intensity 1^ where the nonsequential 
double ionization is expected to be maximum and the intensity where the 
double ionization is complete (see Fig. [I]). In comparison with the helium case, 
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we highlight the common features but also the differences generated by the two 
nuclei, and in particular how the predictions depend on the distance between 
the two nuclei. 



2. Hamiltonian models for H2 



We consider a model for H2 where the two nuclei are fixed [41ll21| . motivated 
by the fact that the time scale of electron motion during the laser pulse is 
short compared to the one characterizing the motion of the heavy nuclei. More 
specifically, we consider the following Hamiltonian model of the H2 molecule (in 



atomic units) with soft Coulomb potentials [41| 
nil ,\ Hp^ II 

1 



IRII 



+ (x + y) • E (t) 



f 



1 



1 



|R/2-x|| 
1 



|R/2 + x|| 
1 



(1) 



f 



|R/2-y|r + l ^||R/2 + y|| 

where x and y are the positions of the two electrons in a <i-dimensional space, 
and p x , p y are their canonically conjugate momenta, and • denotes the Eu- 
clidean scalar product in M. d and ||-|| its corresponding norm. Here we consider 
d = 1, 2 and 3 for the computations, but we mainly consider d = 1 for the phase 
space analysis since we will see below that it reproduces the same characteristic 
features as the higher dimensional cases. The energy is initially fixed at the 
ground state £ g = —1.16 a.u. 
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43j. The linearly polarized laser field is a 



sinusoidal pulse modulated by an envelope, i.e. E(£) = e\Eof(t) smut where 
ei is a unit vector characterizing the direction of the polarization, Eq is the 
maximum amplitude and u> the laser frequency chosen at u> = 0.0584 a.u. which 
corresponds to a wavelength of 780 nm. The pulse envelope f(t) is chosen as 
a trapezoidal function with 2-4-2 laser pulse shape (the ramp-up lasts two cy- 
cles, the plateau four and the ramp-down two). We assume the two nuclei are 
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Figure 1: Double ionization probability for Hamiltonian ([TJ for the d = 1 (dots) ci = 2 (stars) 
and d = 3 (triangles) models as a function of the intensity of the field / for w = 0.0584. The 
vertical lines indicate the laser intensities (continuous green line) ps 3.60 x 10 14 W ■ cm 
and (dashed green line) a 3.60 x 10 14 W • cm where our dynamical analysis predicts the 
maximum of nonsequential double ionization (related to two different approximations), and 
(in blue) the intensity » 7.10 X 10 15 W • cm where the double ionization is expected to 
be complete. 



5 



aligned with the laser field and the distance, denoted R, between them is fixed 
during the whole range of time we will consider [21] , so that R = Rei where 
R = 1.4 a.u. [43 1 . 



Typical ionizing trajectories of Hamiltonian ([T]) show mainly two qualitatively 
different routes to double ionization (see Fig. [2]) : nonsequential double ioniza- 
tion (NSDI), where the two electrons leave the core (inner) region at about the 
same time, and sequential double ionization (SDI), where one electron leaves 
the inner region long after the other one has ionized. 

Using a large assembly of initial conditions, we compute the double ionization 
probability as a function of the laser intensity I (related to Eq by /(W-cm~ 2 ) = 
3.521 x 10 16 £;o(a.u.) 2 ) for d= 1, 2 and 3 (see Fig.[l]). This probability curve has 
the shape of a knee, a hallmark of double ionization. Even if the probability of 
double ionization is much lower for the three or the two dimensional model, one 
can see that characteristic features (the localization of the knee and the intensity 
where double ionization is complete) happen at about the same intensities. We 
restrict ourselves to the one-dimensional case in what follows in order to obtain 
some information on these intensities and the mechanisms behind NSDI and 
SDI. 

We first analyze the dynamics of Hamiltonian ([1]) without the field {Eq = 0). 
This Hamiltonian system has two degrees of freedom, and therefore lends itself 
very well to an analysis using Poincare sections. Here we use a different tool 
to analyze the dynamics, by considering the linear stability properties such as 



obtained by the finite-time Lyapunov (FTL) exponents [4J, |45j . These FTL ex- 
ponents are obtained by integrating the tangent flow together with the equations 
of motion for X = (x, y,p x ,p y ) '■ 

X = F(X,f), (2) 
j = DF(X.,t)J, (3) 

where Eq. ^ are the equations of motion, and Eq. ([3]) is the tangent flow where 
DF(X,t) is the matrix of variations of the generalized velocity field F at the 
point X and time t, i.e. DFij = dFi/dXj. The initial condition for the inte- 
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Figure 2: Example of sequential (upper panel) and nonsequential (lower panel) double ion- 
ization trajectories of Hamiltonian JTJ (for the d = 2 model) for u> = 0.0584 and / = 
3 X 10 15 W ■ cm -2 . The positions along the polarization axis x\ (continuous blue curve) 
and y\ (continuous red curve) of each pair of electrons are plotted versus time, as well as the 
positions in the transverse direction X2 (dashed cyan curve) and yi (dashed magenta curve). 
The pulse shape function f(t) (black curve) and the laser excitation E(t) (dotted black curve) 
are also represented. The amplitude of the shape function and of the laser excitation are not 
representative of the actual conditions. 
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gration of the tangent flow is Jo = I4, the four dimensional identity matrix. 
The (maximum) FTL exponent at time t for the initial conditions Xo is equal 
to Z(i;Xo) = log |A(t;Xo)|/t where A(i;Xo) is the eigenvalue of the Jacobian 
matrix J at time t with the largest norm. The way to analyze the dynamics 
using these exponents is to represent maps of FTL exponents as functions of the 
initial conditions Xo at a fixed time t. These maps (called FTL maps) quantify 
the (linear) instability of some regions and highlight invariant objects. A typical 
FTL map is depicted in Fig. [3] for Hamiltonian |T]) without the field. It clearly 
displays strong and global chaos by showing fine details of the stretching and 
folding of trajectories [44J]. A closer inspection shows some regions associated 
with a very small FTL exponents (blue regions and lower panel of Fig. [3]) . A 
Poincare section reveals the presence of an elliptic island which surrounds an 
elliptic periodic orbit O e , represented on Fig. 3] with a black curve. However 
such regions are too small to play a major role in the dynamics of a typical 
trajectory. In fact, the motion without the field is guided by four weakly hy- 
perbolic periodic orbits. These four important periodic orbits are denoted Xt \, 
O x ,2 1 O y> i and O y .2 1 and the projections of O x ,i and O y ,i are displayed in 
Fig. [U These four periodic orbits are obtained from one single periodic orbit 
through the symmetries of the equations of motion: (x, y,p x ,p y ) >-> (y, 
and (x, y,PxiPy) (-1, —y, —p x , ~Py)- It also means that the representation of 
these orbits on the plane (y,p y ) can be deduced from the one on the plane (x,p x ) 
by inverting the coordinates x and y. The two periodic orbits O x .\ and Xi i 
(stretched in the x-direction) have an outer projection in the (x,p x ) plane and 
an inner one in the (y,p y ) plane (see the red curves in Fig. 2J. By symmetry, 
the same holds for O y ^\ and O a ,2, swapping the role of x and y. Figure [5] repre- 
sents the distance of a typical trajectory with the four periodic orbits mentioned 
above. It shows that at each time, the trajectory is close to one of these orbits 
with rapid transitions between them. This means that the motion on each of 
these periodic orbits, and consequently of a typical trajectory, is composed of 
one electron moving close to the nucleus (the inner electron) and the other one 
further moving away (the outer electron), with quick exchanges of the roles of 
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each electron. 

The reason why these hyperbolic periodic orbits are important is that they 
are short and weakly hyperbolic periodic orbits (their Greene residue [46| is equal 
to 1.1). Hence a typical trajectory passing nearby one of these orbits mimics (by 
continuity) the motion on the periodic orbit and then slowly escapes (depending 
on its escape rate and hence on its stability properties) from this orbit following 
its unstable manifold. 

Single ionization. By switching on the field, the outer electron is picked up 
and swept away from the nucleus. Consequently, its effective Hamiltonian is 
composed of its kinetic energy and the interaction potential of the laser field : 

Hx = ^+E xf{t)smwt (4) 

We notice that the trajectories of Hamiltonian Ti\ can be explicitly computed. 
They are composed of a linear escape from the nuclei (at time to) modulated 



by the action of the field 



22j | (see trajectories after collision in Fig. 



Sequential double ionization (SDI). Since the outer electron is far from the core 
region, the effective Hamiltonian for the inner electron contains the potential of 
the nuclei and the interaction with the laser field : 

P 2 

H 2 = y+y-Bosinwt 

(5) 



^(y-i?/2) 2 + l ^(y + R/2f + 



1 



In the absence of the field (1 = 0), Ti.2 is also integrable and the inner electron is 
confined on a periodic orbit, for almost all the initial conditions. By switching 
on the field, the dynamics associated with Ti.2 is significantly affected. A contour 
plot of y after two laser cycles and a Poincare section are represented on Fig. [6] 
for I = 3 X 10 15 W • cm~ 2 . We notice the following features : First, most of the 
regular trajectories evolve on invariant tori. Second, most of the inner tori are 
broken by the field, creating a bounded chaotic region between the two nuclei. 
These two types of trajectories (regular and bounded chaotic) do not ionize. 
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Figure 3: Upper panel : FTL map of Hamiltonian JTJ without the field at time t = 50 a.u. 
in the plane (x,p x ) with y = 0. The full circles indicate (with the same color code) the 
intersections of the periodic orbits of Fig. [4] with the Poincare section at y = : The 
full diamonds indicate the symmetrical orbits from the one of Fig. [4] with the symmetry 
(x,y,p x ,p y ) i — ► (—x,—y,—p x ,—py). Lower panel : Enlargement of the FTL map around 
a region indicated by a square in the upper panel, and Poincare sections of some regular 
trajectories in the elliptic island. 
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Figure 5: Upper panel : Positions x (red online) and y (blue online) of the two electrons as a 
function of time of a typical trajectory of Hamiltonian JTJ with d = 1 and without the field. 
Lower panel : Distance in phase space of the trajectory depicted in the upper panel versus 
time to the four periodic orbits O x l, O x ,2, O y ^\, O b ,2- For each distance, the color code and 
line style follow the one in Fig. [4] the dashed curves (magenta and cyan) correspond to the 
distance to the symmetrical periodic orbits O x ,2 and O y< 2 respectively. The horizontal dashed 
green lines indicate the positions of the nuclei. Time is in atomic units. 
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Third, the outermost invariant tori are broken since the motion is unbounded 
sufficiently far away from the nuclei. These are the ionizing trajectories. During 
a full laser cycle, we expect two types of non-ionizing trajectories represented 
on Fig. [7] : A typical motion in between the nuclei (upper panel) and another 
surrounding the two nuclei (lower panel). This is what is observed on typical 
trajectories like the one represented on the lower panel of Fig. [5] 

By varying the intensity, the picture of the phase space of Hamiltonian 1(5]) 
evolves in the following way : If the laser intensity I is too small, then the phase 
space is surrounded by invariant tori (with a bounded chaotic region between the 
two nuclei) and no sequential double ionization can occur. The sequential double 
ionization probability depends the on the size of the regular region, and hence 
on I. For large values of the intensity, the regular zone vanishes and after some 
critical intensity, after which there are no invariant tori bounding the motion of 
the inner electron. Consequently, only sequential double ionization is observed 
in this high intensity regime. More specifically, the size of the region where the 
trajectories are bounded is quantified by y rn (defined as the maximum position 
y of the outermost invariant torus of Hamiltonian (|5|) on the Poincare section). 
A numerical approximation to y m is obtained by integrating trajectories at 
p y = (since the domain is approximately symmetric with respect to p y = 0) 
and monitoring the ones which do not keep a bounded value for the position 
after the duration of the laser pulse. This numerical estimate of y m versus the 
intensity / is represented on Fig. [5] (continuous curve). 

A rough approximation to y m — y m {Eq) is given implicitly by the value 
where the potential of Hamiltonian (|5|) is locally maximum, i.e. : 

p Vm - R/2 , Vrn + R/2 
° = 1 2 \W 1 2 \W' { ' 

({y m - R/2) 2 + lj ({y m + R/2) 2 + lj 

independently of the laser frequency. The solution of the above equation is rep- 
resented as a function of the laser intensity / in Fig. [9] (dashed curve). The 
observation is that the theoretical prediction captures the right order of mag- 
nitude, but for u! = 0.0584, this approximation is quantitatively not accurate 
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Figure 6: Upper panel : Contour plot of the electron location y(t) after two laser cycles in the 
plane of initial conditions (yo,P y ,o) of Hamilionian J5j for I = 3 X 10 15 W ■ cm -2 and ui = 
0.0584. The color code is on a logarithmic scale. Lower panel : Poincare section (stroboscopic 
plot with a period of one laser cycle) of some trajectories for the same Hamiltonian as the 
upper panel. The black dots indicate the positions of the nuclei. 
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Figure 7: Two types of trajectories of Hamiltonian J5} for / = 10 15 W • cm -2 and u> = 0.0584 
during a laser cycle : One inside the inner chaotic region (upper panel), and one on an outer 
invariant torus (lower panel). The horizontal dashed lines (green online) indicate the positions 
of the two nuclei. 
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Figure 8: Two typical trajectories of Hamiltonian {TJ for / = 10 15 W • cm -2 and u) = 0.0584 
for initial conditions in the ground state energy of the H2 molecule. The horizontal dashed 
lines (green online) indicate the positions of the two nuclei. 
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Figure 9: Numerical estimates of y m as a function of the laser intensity / for Hamiltonian JTJ 
for w = 0.0584 (continuous blue curve). The qualitative approximation of y m given by Eq. (JBJl 
(dashed curve) and the same qualitative approximation divided by two (thin continuous curve) . 
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since the estimates are about twice the actual values. Of course the quality of 
the approximation depends on the chosen frequency lu. As the frequency to is 
decreased toward zero, the approximation becomes better. 
Equation ^ has a solution for / < 1.09 x 10 16 W • cm~ 2 . For higher intensities, 
we expect complete double ionization. This is the case of Fig. [TJ However, this is 
an upper bound since complete double ionization appears to happen at smaller 
intensities. This emerges from the numerical computation of y m where the size 
of the regular zone drops to zero for intensities larger than 7.10 x 10 15 W-cm~ 2 . 
The origin of this phenomenon is explained by looking at the central periodic 
orbit which organizes the bounded motion. This periodic orbit has the same 
period as the field (and therefore intersects the Poincare section at only one 
point). It is located at y = 0. In order to determine its momentum, we use 



a Newton- Raphson method for the determination of periodic orbits [441 ] . The 
results are represented on Fig. []JjJ We see that for a wide range of intensities 
the momentum p y of the periodic (on the Poincare section) does not change 
significantly. Around I = 7 x 10 15 W • cm -2 it starts decreasing and drops to 
— oo at /(*) = 7.10 x 10 15 W • cm 2 . At this value, complete sequential double 
ionization occurs, in good agreement with Fig. [TJ 

Nonsequential double ionization (NSDI). As noted before, when the field is 
turned on, its action is concentrated on only one electron, the outer one, as 
a first step. The field drives the outer electron away from the nuclei leaving 
the inner electron nearly unaffected by the field (or at least bounded if it is 
between the two nuclei) because its position remains small. From the recollision 



process 



22, 



2l| , the outer electron might come back close to the nuclei during the 
pulse plateau. In this case, it transfers a part of its energy to the inner electron 
through the electron-electron interaction term (a recollision). From then on, 
two outcomes are generically possible : If the energy brought in by the outer 
electron is sufficient for the other electron to escape from the inner region (as in 
the lower panel of Fig. [21) , then it might ionize together with the outer electron. 
The maximum energy E x of the outer electron when it returns to the inner 
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Figure 10: Momentum of the central periodic orbit (on the Poincare section) of Hamilto- 
nian (0 as a function of the laser intensity. The vertical line on the lower panel indicates 
the intensity JW = 7.10 X 10 15 W ■ cm -2 such that for / > /W, complete unhindered SDI is 
expected. The black dots indicate the positions of the two nuclei. Inset : Projection of the 
central periodic orbit at / = 7 X 10 15 W ■ cm~ 2 in the (y, p H )-plane. 



region is obtained from Hamiltonian ((4]) and is £ x = kU v where U p = Eq/(Auj 



is the ponderomotive energy and K = 3.17314 . . . 22j,|47|. We complement the 
recollision scenario (which focuses on the outer electron) by providing the phase 
space picture of the inner electron as accurately described by Hamiltonian ([5]) : 
In order to ionize the inner electron in the most efficient way, the energy brought 
back by the outer electron has to be of order of the energy between the center 
(y = 0) and the boundary of the bounded inner region (y — y m ) of the phase 
space of Hamiltonian H2 (see Fig- [6j) . This energy difference is equal to : 

A£y = 2 l 



1 



^{y m -R/2f + l 



= ■ (7) 

'{y m + R/2f + l 

The equal-sharing relation which links the classical picture of the outer elec- 
tron x with the one of the inner electron y, 



y 2 8uj 2 



(8) 
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defines (through an implicit equation) the expected value of the field Eq and 
hence the corresponding intensity I^ c ' for maximal NSDI because it describes 
the case when each outer electron brings back enough energy to ionize the 
innermost electrons, while keeping enough energy to potentially remain ionized 
itself. However, fulfilling this energy requirement does not guarantee NSDI : The 
outcome depends on the number and the efficiency of recollisions. In order to 
estimate 1^ , we first use the qualitative approximation of y m given by Eq. ([5]) . 
For oj — 0.0584, this intensity is equal to 3.60 x 10 14 W • cm~ 2 which is in good 
agreement with the maximum of the knee in Fig. [1] (green dashed curve) . Now 
we consider the numerically determined values for y m (see Fig. [9j. By solving 
Eq. ©, the intensity 1^ is found at 2.60 x 10 14 W • cm 2 (green continuous 
curve in Fig. [T]) in qualitative agreement with the above value obtained using 
the qualitative determination of y m . 

Next we investigate how this critical intensity 1^ varies with the laser fre- 
quency and highlight how this intensity varies with the distance R between the 
two nuclei. For that purpose, we expand the equal sharing relation {5} using to 
small parameters: R 2 /8 and rjo = 2cj/ v / k (which are assumed to be of the same 
order). The leading order of Eq. ([7]) is 



and the corresponding expansion for the approximation y m — y m (Eo) given by 





(24 + l) 3/2 ' 



If we denote rj = + 1, then Eq. © is expanded into 




An expansion of a solution of the above equation is given by 




which translates into an expansion for the critical amplitude 



E, 



'o 





20 



We notice that the expansion of Eq' at R = is the same as the one found 



Ac) 

for the helium atom [40j since the Hamiltonian model for the helium atom is 
the R = limit of Hamiltonian flj) . As the distance R between the two 
nuclei increases, the location of the knee is displaced toward lower intensities, 
and this displacement is proportional to R 2 . This is actually what is observed 
when we compared the values of 1^ for the helium atom 40 1 (7"( c )(i? = 0) = 
4.60 x 10 14 W • cm" 2 for uj = 0.0584) and for the H 2 molecule. For u = 0.0584, 
the approximate value of 1^ given by Eq. (J9]) for R — 1.4 is 3.40 x 10 14 W-cm~ 2 
which is in good agreement the numerical solution of Eq. (|8|) . 



3. Conclusion 

The classical picture for the sequential and nonsequential double ionization 
of H2 is obtained by complementing the recollision scenario for the outer electron 
with the phase space picture of the inner electron. By finding the organizing 
principles of the classical dynamics, we arrived at two predictions for the char- 
acteristic intensities of the knee-shaped double ionization probability curve : 
the intensity after which the complete double ionization is expected, and the in- 
tensity where the non-sequential double ionization is predicted to be maximum. 
Very good agreement is found by comparing our predictions with the direct inte- 
gration of a large assembly of trajectories. It should be noted that this scenario 
follows closely the one of the helium atom, implying the existence of a general 
mechanism for nonsequential double ionization of atoms and molecules. 

Based on these findings, we believe that advanced methods of classical me- 
chanics and the diagnostics we use here are worth pursuing further because 
classical mechanics, with its advantageous scaling with system size, may well 
become a useful tool for understanding some aspects of molecular systems too 
complex for a full quantal treatment with contemporary computing resources. 
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